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Abstract 

Two-point functions, mean-squared fluctuations, and tlie vacuum expectation value of the 
energy-momentum tensor operator are investigated for a massive scalar field with an arbitrary 
(~| , curvature coupling parameter, subject to a spherical boundary in the background of de Sitter 

spacetime. The field is prepared in the Bunch-Davies vacuum state and is constrained to 
satisfy Robin boundary conditions on the sphere. Both the interior and exterior regions 
are considered. For the calculation in the interior region, a mode-summation method is 
employed, supplemented with a variant of the generalized Abel-Plana formula. This allows 
us to explicitly extract the contributions to the expectation values which come from de Sitter 
spacetime without boundaries. We show that the vacuum energy-momentum tensor is non- 
diagonal with the off-diagonal component corresponding to the energy flux along the radial 
direction. With dependence on the boundary condition and the mass of the field, this flux 
can be either positive or negative. Several limiting cases of interest are then studied. In terms 
of the curvature coupling parameter and the mass of the field, two very different regimes 
are realized, which exhibit monotonic and oscillatory behavior of the vacuum expectation 
values, respectively, far from the sphere. The decay of the boundary induced expectation 
values at large distances from the sphere is shown to be power-law (monotonic or oscillating) , 
independent of the value of the field mass. The expressions for the Casimir densities in the 
exterior region are generalized for a more general class of spherically-symmetric spacetimes 
inside the sphere. 



PACS numbers: 04.62. -Fv, 03.70. -Fk, ll.lO.Kk 



1 Introduction 

De Sitter (dS) spacetime is one of the simplest and most interesting spacetimes allowed by 
general relativity. Quantum field theory in this background has been extensively studied during 
the past two decades. Much of the early interest was motivated by the questions related to 
the quantization of fields on curved backgrounds. dS spacetime has a high degree of symmetry 
and numerous physical problems are exactly solvable on this background. The importance of 
this theoretical work increased by the appearance of the inflationary cosmology scenario fT]. 
In most inflationary models, an approximately dS spacetime is employed to solve a number of 
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problems in standard cosmology. During an inflationary epoch, quantum fluctuations in the 
inflaton field introduce inhomogeneities which play a central role in the generation of cosmic 
structures from inflation. More recently, astronomical observations of high redshift supernovae, 
galaxy clusters, and cosmic microwave background [2j indicate that at the present epoch the 
universe is accelerating and can be well approximated by a world with a positive cosmological 
constant. If the universe would accelerate indefinitely, the standard cosmology would lead to an 
asymptotic dS universe. Hence, the investigation of physical effects in dS spacetime is important 
for understanding both the early universe and its future. Another motivation for investigations 
of dS-based quantum theories is related to the holographic duality between quantum gravity 
on dS spacetime and a quantum field theory living on a boundary identified with the timelike 
infinity of dS spacetime |3]. 

In dS spacetime, the interaction of fluctuating quantum fields with the background grav- 
itational field gives rise to vacuum polarization. The Casimir effect presents another type of 
vacuum polarization induced by the presence of boundaries. This effect is among the most 
striking macroscopic manifestations of non-trivial properties of the quantum vacuum. It has 
important implications on all scales, from subnuclear to cosmological. The reflecting boundaries 
alter the zero-point modes of a quantized field and shift the vacuum expectation values of quan- 
tities, such as the energy density and stresses. As a result, forces arise acting on constraining 
boundaries. The particular features of these forces depend on the nature of the quantum field, 
the type of spacetime manifold, the boundary geometry, and the specific boundary conditions 
imposed on the field. Since the original work by Casimir many theoretical and experimental 
works have been done on this problem (see, e.g., Ref. ^ and references therein). 

In the present paper we consider a problem with both types of vacuum polarization. Namely, 
we evaluate the vacuum expectation values for the field squared and the energy-momentum ten- 
sor of a scalar field with general curvature coupling parameter induced by a spherical boundary 
on the background of {D + l)-dimensional dS spacetime. Historically, the investigation of the 
Casimir effect for a spherical shell was motivated by the Casimir semiclassical model of an elec- 
tron. In this model Casimir suggested [5j that Poincare stress, to stabilize the charged particle, 
could arise from vacuum quantum fluctuations and the flne structure constant could be deter- 
mined by a balance between the Casimir force (assumed attractive) and the Coulomb repulsion. 
However, as has been shown by Boyer the Casimir energy for a perfectly conducting sphere 
is positive, implying a repulsive force. This result was later reconsidered by a number of authors 
[7j. More recently new methods have been developed for the investigation of the Casimir effect 
including direct mode summation techniques and the zeta function regularization scheme [8], 
semiclassical methods [9] in the framework of the Gutzwiller trace formula, the optical approach 
|10j , worldline numerics [IT] , the path integral approach [12j , methods based on scattering theory 
|13] . numerical methods based on evaluation of the stress tensor via the fluctuation-dissipation 
theorem |14j (for a review see Refs. [U [15]). The Casimir effect for a spherical shell in an ar- 
bitrary number of dimensions is analyzed in Refs. [16\ ITT] for a massless scalar field satisfying 
Dirichlet and a special type of Robin (corresponding to the electromagnetic TM modes) bound- 
ary conditions using the Green's function method and in Refs. |18l [19] for the electromagnetic 
field and massless scalar and spinor fields with various boundary conditions on the basis of the 
zeta regularization technique. The case for a massive vector field has been discussed in Ref. |20j . 

The investigation of the energy distribution inside a perfectly reflecting spherical shell was 
made in Ref. [21]. The distribution of the other components for the energy-momentum tensor 
of the electromagnetic field inside and outside a shell and in the region between two concentric 
spherical shells is studied in Refs. [22l [23] (see also Ref. [23]). The investigation of the elec- 
tromagnetic energy density near a conducting boundary was first carried out by DeWitt |25| . 
The vacuum expectation values for the energy-momentum tensor of a massive scalar field with 
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general curvature coupling parameter and obeying the Robin boundary condition on spheri- 
cally symmetric boundaries in (D + l)-dimensional spacetime are investigated in Ref. [26j. The 
Casimir densities for spherical boundaries in the background of global monopole and Rindler-like 
spacetimes have been discussed in Refs. [27J and [28j respectively. 

Previously, the Casimir stresses for spherical boundaries on the background of dS spacetime 
have been investigated in Ref. [29] for a conformally coupled massless scalar field. In this last 
case the problem is conformally related to the corresponding problem in Minkowski spacetime 
and the vacuum characteristics are generated from those for the Minkowski counterpart, just 
by multiplying with the conformal factor. As it has been shown in Refs. [30] for the geometry 
of flat boundaries, qualitatively new features arise in the case of non-conformally coupled fields 
(for the case of a minimally coupled massless field see Ref. |31j). (For flat backgrounds, see 
also Ref. |32j.) The curvature of the background spacetime decisively influences the behavior 
of boundary-induced vacuum expectation values at distances larger than the curvature scale. 
Recently, the topological Casimir effect in dS spacetime with toroidally compactified spatial 
dimensions has been investigated in Ref. [33j • 

We have organized the paper as follows. In the next section the Wightman function is 
evaluated inside a spherical boundary in dS spacetime for a scalar field with general curva- 
ture coupling parameter and with Robin boundary condition on the sphere. Among the most 
important quantities describing the local properties of a quantum field and the corresponding 
quantum back-reaction effects are the expectation values of the field squared and of the energy- 
momentum tensor. These quantities in the interior of a sphere will be investigated in Sees. [3] 
and m The Wightman function for the region outside a spherical boundary is considered in 
Sec. [5l The vacuum expectation values of the field squared and the energy-momentum tensor 
in this region are discussed in Sec. [H Sec. [7] generalizes the results for the exterior region to a 
more general class of spherically-symmetric spacetimes inside a spherical boundary. Section [8] 
contains a summary of the work. In Appendix [X] the expression for the Wightman function in 
boundary-free dS spacetime is derived by making use the corresponding mode sum. In Appendix 
[B| for the Wightman function, we explicitly demonstrate the limit to the geometry of a spherical 
boundary in Minkowski spacetime. Appendix O sketches the corresponding calculation of the 
Green's function for this problem. 



2 Wightman function inside a sphere 

We consider a quantum scalar field (p{x) on a (D + l)-dimensional dS spacetime background 
described in inflationary coordinates. The latter are most appropriate for cosmological applica- 
tions. The spatial part of the line element we will write in terms of the hyperspherical coordinates 
(r,i?,</>) = {r,ei,e2,...en,(l)),n = D-2: 

ds^ = dt^ - e2*/°((ir2 + r2(if]2^_i), (2.1) 

where d^^j^_-^ is the line element on a (-D — l)-dimensional sphere with unit radius. In what 
follows, in addition to the synchronous time coordinate, t, we will also use the conformal time, 
r, defined as r = — ae~*/", — oo < r < 0. In terms of this coordinate the line element takes the 
conformally flat form: 

ds^ = a^T-"^ [dr^ - dr^ - r^dfll^^) . (2.2) 

Note that the parameter a is related to the cosmological constant A through the expression 
a = D{D -1)/{2A). 

The dynamics of a massive scalar field are governed by the equation [34] 

{ViV + + (R)ip = 0, (2.3) 
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where is the covariant derivative operator, R = D{D + l)/a^ is the Ricci scalar for dS 
spacetime, and ^ is the curvature coupling parameter. The special values of this parameter 
^ = and ^ = = {D — 1)/4D correspond to minimally and to conformally coupled fields. 
The importance of these two special cases comes from the fact that, in the massless limit, 
the corresponding fields mimic the behavior of gravitons and photons, respectively. Our main 
interest in this paper is the study of the changes in the vacuum expectation values (VEVs) of 
the field squared and the energy-momentum tensor induced by a spherical shell with radius a 
centered at the origin in dS spacetime. We assume that on the sphere the field obeys Robin 
boundary condition 

(A + Bdr)ip{x) =0, r = a, (2.4) 

with constant coefficients A and B, in general, different for the inner and outer regions. The 
results for Dirichlet and Neumann boundary conditions are obtained as special cases. Robin 
boundary conditions are an extension of the ones imposed on perfectly conducting boundaries 
and may, in some geometries, be useful for depicting the finite penetration of the field into the 
boundary with the "skin-depth" parameter related to the Robin coefficient )35j . These types of 
conditions naturally arise for scalar and fermion bulk fields in braneworld models. 

As the first step in the investigation of the VEVs we will evaluate the Wightman function 
W{x, x') = {0\ip{x)(p(x')\0) , where |0) stands for the vacuum state (for the Wightman function in 
the geometry of spherical boundaries in the background of a constant negative curvature space 
see Ref. [36]). In order to do that we employ the mode sum formula 

W{x,x') = Y,M^Mi^'), (2.5) 

with {ipo-{x),ip*{x)} being a complete set of solutions to the classical field equation, specified 
by a set of quantum numbers a, satisfying the boundary condition (j2.4p . In accordance with 
the spherical symmetry of the problem under consideration, the angular dependence of the 
mode functions is given by the spherical harmonic of degree / (see Ref. [37]), Y{mp] where 
nip = (mo = /, mi, . . . , m„), / = 0, 1, 2, . . ., and mi, m2, . . . , m„ are integers such that 

^ m„_i ^ m„_2 ^ • • • ^ mi ^ /, 

-rrin-i ^ m.„ ^ m„_.i. (2.6) 

Presenting the mode functions in the form (paix) = T {T)TZ{r)Y (mp; i!) , cf)) , from the field 
equation (j2.3p it follows that the time and the radial coordinate dependences are given in terms 
of cylinder functions as: 

= ri-^/'[6iJ^(Ar) + 62y^(Ar)], 
r(r) = r,""^' c,i/^)(Ar?), ^ = \t\, (2.7) 

where Ji_i{z) and Y^{z) are the Bessel and Neumann functions respectively, Hu\z), j = 1, 2, are 
the Hankel functions (we use the notations from Ref. [38] ) . The orders of the cylinder functions 
in Eq. ()2.7p are defined as: 

^ = / + L>/2-l, 

u = [D^/4-D{D + l)C-m'^a'^Y^^ . (2.8) 

Note that i' is either real and nonnegative or purely imaginary. For a conformally coupled 
massless field u = 1/2 and the Hankel functions in Eq. (j2.7p are expressed in terms of elementary 
functions. 
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Different choices of the coefficients cj in the expression for the function T{t) correspond 
to different choices of the vacuum state in dS spacetime. The choice of the vacuum state is 
among the most important steps in construction of a quantum field theory in a fixed classical 
gravitational background. dS spacetime is a maximally symmetric space and it is natural to 
choose a vacuum state having the same symmetry. In fact, there exists a one-parameter family 
of maximally symmetric quantum states (see, for instance, Ref. [39] and references therein). 
Here we will assume that the field is prepared in the dS-invariant Bunch-Davies vacuum state 
|40) for which C2 = 0. Among the set of dS-invariant quantum states the Bunch-Davies vacuum 
is the only one for which the ultraviolet behavior of the two-point functions is the same as in 
Minkowski spacetime. 

First we consider the region inside the spherical shell. From the regularity condition at the 
origin it follows that for this region 62 = and the mode functions realizing the Bunch-Davies 
vacuum state are written in the form 

'^-(^) = C„^^^Hi^\>^v)M>^r)Y{mp;i),cj)). (2.9) 

From the boundary condition ()2.4p . it follows that the eigenvalues for A have to be solutions to 
the equation 

AJ^,{Xa) + B\aJ'^{\a) = 0, (2.10) 

where the prime means the derivative with respect to the argument of the function. In Eq. (j2.10p 
and in what follows we use the notations 

A = A + {1- D/2)B/a, B = B/a. (2.11) 

For real A, B and /i > —1, all roots of Eq. (I2.10j) are simple and real, except the case A/B < 
when there are two purely imaginary zeros (see, e.g., Ref. [H]). We will denote by 2; = A^^fc, k = 
1,2,..., the zeros of the function AJ^{z) + BzJ'^{z) in the right half-plane, assuming that they 
are arranged in ascending order. So, for the eigenvalues of A one has A = X^^k/o,- Now we see 
that the set of quantum numbers a is specified to cr = (A;, /, mi, . . . , m„). Note that for Dirichlet 
boundary condition A = \, B = and for Neumann boundary condition A = [\ — D/2)B. 
The coefficient in Eq. (j2.9p is determined from the orthonormalization condition 



; j d^X.Md'^aix)^ r^l^ix) = 6^^,, (2.12) 



where the integration goes over the region inside the sphere and Jo-o-' is the Kronecker delta. 
Substituting the functions from Eq. (j2.9p . the integration over the angular variables is performed 
by making use of the normalization integral for the spherical harmonics: 



/ 



\Y{mk;^A)fdn = N{mk), (2.13) 



where the explicit form for N{mk) can be found in Ref. ^7J. By using also the relation (for u 
purely real or imaginary) 

Hi'\Xv)H;^J'{Xr^) - i/i?(A7?)i/p(Ar/) = -iLe-(--*)-/2, (2.14) 

nXr] 

for the normalization coefficient we get 

, _ ^AT^(Aa)e-(--*)-/^ 
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with the notation 

2; 

'^'^"^ = (z2-^2)j2(,) + ,2j;2(,)- (2-16) 

Substituting the eigenfunctions into the mode sum formula ()2.5p , for the Wightman function 
one gets {n = D — 2) 

00 

xY,zT^{^)H^u\^^lo^)H'${zr^' /a)J^{zr/a)J^{zr' /a)\^^^^^. (2.17) 
fc=i 

Here Cp{x) is the Gegenbauer or ultraspherical polynomial of degree p and order g, 5^) = 
2-K^/'^/T{D/2) is the surface area of a unit sphere in D-dimensional space, and is the angle 
between the directions {'&,4>) and (i?', (/>'). In deriving Eq. (j2.17p . we have used the addition 
theorem [37] for the spherical harmonics: 

^Yir^ = ^C^lcosS), (2^18) 

nip 

where the sum is taken over the integer values m^, p = 1, 2, . . . , n, in accordance with Eq. (j2.6p . 

Before proceeding to the evaluation of the Wightman function, we comment about the realiz- 
ability of the Bunch-Davies vacuum state. It is well known that in dS spacetime without bound- 
aries the Bunch-Davies vacuum state is not a physically realizable state for Rei^ ^ D/2. The 
corresponding Wightman function contains infrared divergences arising from long-wavelength 
modes. As it has been shown in Ref. [32], these divergences lead to inconsistencies with Einstein 
equations and they cannot arise through dynamical evolution from a state which is initially free 
of such divergences. In the presence of boundaries, the boundary condition imposed on the 
quantized field may exclude these modes and the Bunch-Davies vacuum becomes a realizable 
state. An example is provided by a spherical boundary described above. In the region inside the 
spheres and for boundary conditions with A 7^ 0, there is a minimum value for A, A ^ A„/2,i/<^) 
and the two-point function (|2.17p contains no infrared divergences. 

The eigenvalues X^^k are given implicitly and the Wightman function in the form (I2.17P is 
not convenient for the further evaluation of the VEVs. In order to sum over these eigenvalues 
we use the summation formula [43j 



00 

2VT^(A^,fc)/(A^,fc) = / /(x)dx + -Res,=o/(^) 

tr, Jo ^ 



k=l 

- - r dx^4^ fe-'^"7(xe"*/2) ^ e'^"7(xe-"^/2)l ^(2.19) 

where f{z) is an analytic function on the right half-plane and, for a given function F{z), the 
barred notation is defined as: 

F{z) = AF{z) + BzF'{z). (2.20) 

Formula (j2.19p can be generalized for the case of the existence of purely imaginary zeros of the 
function Jp{z) by adding the corresponding residue term and taking the principal value of the 
integral on the right (see Ref. [l3]). In what follows we assume values ol A/B for which all roots 
A/^,A; are real. 
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As a function f{z) in Eq. (|2.19|) we take 



fiz) = zHi^\zr]/a)Hl^J {ZT]' /a)J^{zr/a)J^izr' /a). (2.21) 

By making use of the properties of the cyhnder functions, after the apphcation of Eq. ()2.19p . 
the Wightman function is presented in the decomposed form: 

Wix,x') = Wdsix,x') + Wi,{x,x'). (2.22) 

Here, the first term in the right-hand side corresponds to the first integral on the right of Eq. 

(EH, 

^ ^ I 1=0 

POD 

X / dXXHi^\Xr])HlV{Xr]')J^{Xr)J^{Xr'), (2.23) 

and 







W,ix,x') = -^|^E(2^ + n)Cr/^(cos^)y^ dz z 



X^^^I^{zr)I^{zr') [l,{zv')K,{zr]) + I^,{zv)K,{zr,')] . (2.24) 
I^{za) 

The term ()2.23p does not depend on the radius of the sphere whereas the second term on 
the right-hand side vanishes in the hmit a — )• oo. From here it fohows that Wdsix,x') is the 
Wightman function for a scalar field in boundary-free dS spacetime. In Appendix |A] we show 
this by direct evaluation. The second term in the right-hand side of Eq. (I2.22P is induced by the 
presence of the spherical shell. In the limit a — t- oo with fixed t, the line element (j2.ip goes to 
the Minkowskian line element in spherical coordinates. In Appendix [B] it is shown that in this 
limit, from Eq. (j2.24p . the Wightman function is obtained for a scalar field inside a spherical 
boundary in Minkowski spacetime. In Appendix C we derive the corresponding causal Green's 
function. 

In a similar way we can evaluate the Wightman function in a general state described by 
the modes (12. 7p without the specification of the coefficients in the linear combination of the 
Hankel functions. After the normalization of the corresponding mode functions, we obtain a set 
of quantum states determined by the ratio c^jcx. In general, the latter may depend on / and A. 
After the application of the generalized Abel-Plana formula, the boundary-induced part in the 
corresponding Wigthman function is presented in a form similar to Eq. (12.24p where now, instead 
of the combination of functions in the square brackets, a more general bilinear combination of 
the modified Bessel functions appears. 



3 VEV of the field squared inside a sphere 

Given the Wightman function, we can proceed to the evaluation of the VEV of the field squared. 
This VEV is among the most important quantities in discussing the phase transitions in the 
early universe and the generation of the cosmic structures from inflation. The VEV of the field 
squared is obtained taking the coincidence limit of the arguments. In this limit the Wightman 
function is divergent and some renormalization procedure is needed. The important point here 
is that for points away from the sphere the divergences are the same as those for dS spacetime 
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without boundaries. As we have akeady extracted the part PVdsC^^; 2;'), the renormahzation is 
reduced to the renormahzation of the VEV in the boundary-free dS spacetime, which is aheady 
done in the literature. In this way, the renormalized VEV of the field squared inside a sphere is 
presented in the decomposed form 

iv') = (9^')dS + {^\, (3.1) 

where (</5^)dS is the VEV in the boundary-free dS spacetime and the term {ip'^)h is induced by 
the sphere. Due to the maximal symmetry of the Bunch-Davies vacuum state the VEV (v'^)ds 
does not depend on the spacetime point. 

The boundary-induced part is directly obtained from the corresponding part in the Wight- 
man function, given by Eq. (|2.24p . By taking into account that 

Cf'in ^ (3,2) 

we get 

(^')b = - \"D-2 E^^ / dxxyij^lf,{rx)K,{xv)[I.{xv) + I-u{xv)]. (3.3) 
In Eq. ([O]) . 

Di = {2l + D- 2) (3.4) 

' ^ ' r{D-i)n ^ ' 

is the degeneracy of the angular mode with given /. The integral representation ()3.3p is valid for 
Rei^ < 1. For large values of x, the integrand in Eq. (13. 3p behaves as g-C'^-'')^. The presence of 
the spherical boundary breaks the dS-invariance and the mean-squared fluctuation of the field 
depends on time. As is seen from Eq. (|3.3p . this dependence appears through the ratios a/rj 
and r/r]. The latter property is a consequence of the maximal symmetry of the Bunch-Davies 
vacuum in the absence of the sphere. Note that a/r] and r/r] are the proper radius of the sphere 
and the proper distance from the sphere center measured in units of the dS curvature scale a. 
The influence of the gravitational field on boundary-induced quantum effects appears through 
the function Ki,{z) \Iu{z) + I-u{zy\. The latter is a monotonically decreasing function of z for 

^ and exhibits an oscillatory behavior for < and ^ As it will be seen below, this 
feature leads to interesting physical consequences. 

For a conformally coupled massless scalar field, ^ = (,d, m = 0, one has z/ = 1/2 and 

Ki/2{x) [h/2{x) + I-i/2{x)] = l/x. (3.5) 
In this case, from Eq. (j3.3p . we find 

(99% = (r?/a)^-i(^2)M, (3.6) 

where 

W M = ^Z^^U dx If,[rx), (3.7 

is the corresponding VEV inside a spherical boundary in Minkowski spacetime [26]. Of course, 
this result could be obtained directly by using the conformal relation between the problems in 
dS and Minkowski spacetimes. 
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Now let us consider the asymptotics of the boundary-induced VEV. Near the center of the 
sphere the dominant contribution comes from the mode with the lowest orbital momentum / = 
and in the leading order we find 



2. (2a)^-^ r, D-iKD/2-iiax/r,) 



^^"b « - D/2+iVm/9^ / '^^^ r [Iu{x)+I-u{x)] . (3.8) 

Tr^/2+iY{D/2) Jo lD/2-i{ax/V) 

The behavior of the field squared for large values of the sphere proper radius, a/i] » 1, assuming 
that r/rj is fixed, is obtained from Eq. (j3.8p expanding the function Ky{x) [Iu{x) + I-u{x)\ for 
small values of the argument. For x ^ 1 in the leading order one has 



K^{x) [Iy{x) + I-u{x)] w o-^Re 



r(l - v)x'^'' 



(3.9) 



with the notation 



where cr^ = 1 for positive v and cji^ = 2 for imaginary v. The case v = Q should be considered 
separately. The behavior of the VEV is qualitatively different for these two cases. For positive 
one gets 

("P )b - - r{D/2) ^^""^^ (^-^"^ 

r(i - 1") Jo Id/2-i{x) 

In this case the boundary-induced VEV monotonically decreases with increasing proper radius 
of the sphere. For imaginary values of z/, the leading term is in the form 

^ -^!^£r-f^^<^) ln(a/,7) + <^(^)], (3.12) 

where c{u) and (/'{i') are defined by the relation b{u) = c{u)e'^'^^'^\ In this case the VEV exhibits 
a damped oscillatory behavior. 

The VEV of the field squared diverges on the boundary. Surface divergences are well known 
in quantum field theory wtih boundaries and they have been investigated for various geometries. 
Near the sphere the dominant contribution in Eq. (|3.3p comes from large values of x and I. By 
taking into account that for large x and for fixed one has K,^{x) [I^ix] + I-^{x)] ^ 1/x, we 
conclude that near the boundary, to the leading order, one has (</J^)b ~ ijl/oi)^^^ {h^'^)m- By 
using the asymptotic expression for (93^)m, we get 

2 ^ {n/ar-'T{{D -l)/2) 

)b - -(4vr)(^+i)/2(a-r)^-i ^^'^^^ 

where 

KB = 26oB - 1. (3.14) 

In deriving Eq. (|3.13|) we have assumed that (a — r) <^ r]/{ma) which corresponds to small 
proper distances from the sphere compared with the Compton wavelength of the scalar particle. 
As the boundary-free part in the VEV is constant everywhere, we conclude that near the sphere 
the total VEV is dominated by the boundary- induced part. The factor in Eq. (j3.13p indicates 
that the sign of the surface divergence is reversed for Dirichlet as opposed to any other Robin 
boundary condition. 
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4 Vacuum energy- momentum tensor inside a sphere 



Another important quantity characterizing the properties of the quantum vacuum is the VEV of 
the energy- momentum tensor. In addition to describing the physical structure of the quantum 
field at a given point, the energy-momentum tensor acts as the source of gravity in the Einstein 
equations. It therefore plays an important role in modelling a self-consistent dynamics involv- 
ing the gravitational field. (For the renormalization of the Einstein equations in the presence 
of energy-momentum tensor divergences, see Ref. [44]; the demonstration that the finite and 
divergent parts of the Casimir energy obey the equivalence principle appears in Refs. [45 ^ 146 j.) 
Similar to the field squared, the VEV of the energy- momentum tensor is decomposed as 



{Tik) = (Tik) ds + {Tik)h, 



(4.1) 



where the boundary-free part is presented in the form {Tik)dS =const-gik- The latter property is 
a direct consequence of the maximal symmetry of the Bunch-Davies vacuum state. In particular, 
for = 3 dS spacetime the renormalized boundary-free part is given by the expression [40^ |47] 
(see also Ref. [3l]) 



{Tik)ds 



9ik 



327r2a4 



{m^a^ 6^ - l) (3/2 + i^) + (3/2 



In ( 



(6^ - ly + 1/30 + (2/3 - e^m^Q^}, 



(4.2) 



where 'ip{x) is the logarithmic derivative of the gamma- function. For ma ^ 1, to the leading 
order, {Tik)ds ^ C gik / {38Att'^ m'^ a^) , where the coefficient C depends on the curvature coupling 
parameter only. For minimally and conformally coupled fields we have C = 7 and C = —1/5 
respectively. 

For the evaluation of the boundary-induced part in the VEV of the energy-momentum tensor, 
we use the formula 



ik/h 



lim did'i,W,,{x, x') + (e - 1/4) QikViV' - eViVfc - ^Rik {ip\, 



(4.3) 



where Rik = Dgik/a^ is the Ricci tensor for dS spacetime. By taking into account Eqs. (j2.24p 
and (13. 3p . after long but straightforward calculations, the VEVs for the diagonal components 
are presented in the form (no summation over i) 



^2 D j^Q: 

'27r5Da^+ 



(=0 



dx X 



I^{ax) 



{G][I^{rx)]F,{xri) + 2ll{rx)Fl(xr,)] , (4.4) 



where we have introduced the notation 

Other functions in Eq. ()4.4p are defined in accordance with (no summation over k) 



(4.5) 



i- 



D dy 

y 



1 + 



e(2-i?) + 



D-l 



Fuiy), 



V 9 

y y^ . 



F,iy), 



(4.6) 
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for k = 1,2, . . . , D, and 

GO[/(y)] = (1-40 

G\ ifiy)] 

Gl [/(?/)] 
for k = 2, . . . , D, where 



D 



y 



-f{y)f'{y) + {i + ^]ny) 



f'\y) + -f{y)f{y)- i + 
y 



(4e-l)r(y)-^/(y)/'(y) + 

y 



{D/2 - 1)^1 + ^2 



4e-i + 



f{y), 



(z?/2-i)6+^2(i+a) 



(4.7) 



{D - l)y2 



f\y), 



^1 = 4(D - i)e - + 2. 

In addition, the boundary-induced VEV has also nonzero off-diagonal component 

'0 



(4.8) 



D 



1=0 



X [2rxl'^{rx) - {D - 2)1 ^{rx)] [(^ - 1/4) r?a^ + F,{xr^), 



(4.9) 



which corresponds to the energy flux along the radial direction (see below). Note that in the 
formulas above, the components of the energy-momentum tensor are written in the coordinates 
(r, r, -i?, (/)). The components in the system {t,r,'d,(f)) with the synchronous time coordinate will 
be denoted by (7'(^)j)b- We have the relations (no summation over i) 



{T, 



{s)i/ 



(4.10) 



As an additional check of the expressions for the VEV of the energy-momentum tensor, it 
can be seen that the boundary-induced parts fulfill the trace relation 



(r/)b = \d{C - CdWiV + m^] {ip\. 



(4.11) 



In particular, the boundary- induced part in the VEV of the energy-momentum tensor is traceless 
for a conformally coupled massless scalar field. The trace anomaly is contained in the boundary- 
free part only. In addition, the VEVs obey the covariant conservation equation Vk{T^)h = 0. 
For the geometry of interest, this equation is reduced to two equations: 



dv - ^) (7^0 )b - drm\ + i [(Ti)b + {D- l)(r|)b] 



D + l 
V 



i/b 



D - 1 



For a conformally coupled massless field one has = 1/2 and Fu{z) 
it is easily seen that the off-diagonal component vanishes and FQ{y) = - 
k = 1,2, . . . , D. As a result we find (no summation over i) 



,D-1 



(4.12) 

With this 
, F^{y) = for 

(4.13) 



where {T-)j^ is the corresponding VEV for a sphere in Minkowski spacetime |26j . Again, the 
result (|4.13p directly follows from the conformal relation between the problems in dS and 
Minkowski spacetimes. The electromagnetic field is conformally invariant in 4-dimensional 
spacetimes. Hence, the electromagnetic Casimir densities for a conducting spherical shell on 
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background of D = 3 dS spacetime are obtained from those in Minkowski spacetime by using 
Eq. 

Let us denote by Ey^ the boundary-induced part of the vacuum energy in the spatial volume 
V with a boundary dV: 

4^) = j^d^x^{Tl^,%, (4.14) 

where 7 is the determinant of the spatial metric tensor 7^5 with gi^ = (1, —7^,5) and with 
the Greek indices running over 1,2, Now, from the covariant conservation equation 

Vfc(r(^)j)b = with i = 0, it follows that 

= - d^-^x Vhnp{T^^^^)^ + i ^ (i^xV7(r(f )^)b, (4.15) 

where n^, ■^^^npn^ = 1, is the external normal to the boundary dV and h is the determinant 
of the induced metric hp^ = 7^5 — npn^. The first term in the right-hand side of Eq. (j4.15p 
describes the energy flux through the boundary dV . In particular, for a spherical boundary 
with radius rg < a and with the center at r = one has = {a/'r])5^ and the flux term is 
given by —{aro/r])^^^ SD{a/r]){T^^^f^)i^. Now, by taking into account that the proper surface 
area of the sphere with radius ro is given by (aro/r/)^~^5'D, we conclude that the quantity 
{Tq)\, = (a/f/)(7'(-s')o)b) given by Eq. (j4.9p . is the energy flux density per unit proper surface 
area. 

By using the expression (j4.4p for the energy density, for the boundary-induced vacuum energy 
inside the sphere with radius ro one gets 

2ia ^ 'J„ /„{ai/r„) 



D — 2 



{{l-AOF,{xr]/ro)I^{x) 
-F^ixrj/ro) [x'l';^{x) - (x^ + f,^)lf,{x)] }. (4.16) 



In deriving Eq. (j4.16p we have used the formula for the integral involving the square of the 
modified Bessel function [IS]. Related to the surface divergences in the vacuum energy density 



(see below), the integrated energy E^^^^ diverges in the limit ro — )• a. In order to obtain a 
finite result for the vacuum energy inside the sphere, an additional renormalization procedure 
is needed. Our main interest in the present paper are local characteristics of the vacuum away 
from the boundary. As we have noted before, the renormalization of the latter is reduced to 
that for the boundary-free dS geometry. 

The general expressions (j4.4p and (j4.9p are simplified at the sphere's center and near the 
sphere. Near the center, the dominant contribution to the boundary-induced VEVs comes from 
the modes with I = and 1 = 1. To the leading order, for the diagonal components we get (no 
summation over i) 

^ b ~ d72TTF77T7^ / dxxl^- — —G x 



Kr)/ri__T, (ax/ri) r . • 1 1 

+ r V GI,^F,{x) + 2F:{x) , (4.17) 

lD/2-i{ax/i]) L ^ ^ J J 
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where we have defined 



'-'{o) - ^ - '-^(0) 

= Z^i-ii, Gf,)=4(D-l)^-Z) + 2, (4.18) 

for k = 1,2,... ,D. The first and second terms in the figure braces of Eq. (j4.17p come from the 
modes I = 1 and / = 0, respectively. As we could expect, the stresses are isotropic at the center. 
For the off-diagonal component of the VEV one has the following leading term: 



2i-^a-^-i 



^ vr^/2+ir(D/2 + l)yo [lD/2-i{ax)^ Io,2{ax) ^ 

x[{i-l/A)r,d^ + i]F,{xr^). (4.19) 

The first and second terms in the square brackets are the contributions of the modes / = 
and 1 = 1. The off-diagonal component vanishes at the center. The behavior of the energy- 
momentum tensor components for large values of the sphere proper radius, a/i] ^ 1, assuming 
fixed values of r/ry, is obtained from Eqs. ()4.17p and (j4.19p by expanding the function Fy{z) for 
small values of the argument. By using Eq. (j3.9p . similar to the case of the field squared, it 
can be seen that for positive values of v the VEVs decay monotonically like {a/r])^'^~^~'^ . For 
imaginary v, the behavior of the VEVs is damped oscillatory with the amplitude decaying as 
(a/r?)"^-2. 

For points near the surface of the sphere, the dominant contribution to the boundary-induced 
VEVs comes from large values of / and x. By using the uniform asymptotic expansions for the 
modified Bessel functions, it can be seen that the leading terms in the diagonal components for 
a scalar field with non-conformal coupling (^ 7^ ^/j) are related to the corresponding terms for a 
spherical boundary in Minkowski spacetime by (no summation over i) {T-)i^ ~ (r//a)^+^(T?)M- 
These leading terms are given by the expressions (no summation over A;) 



roo 



KD/2-i{ax) Kr)/2{ax) 



DT{iD + l)/2)i^-^D)^B 
2Dn(D+i)/2[a{a-r)/r]] 



^^i/b - 2^7r(^+i)/2a(a-r)^(a/r/)^+i ' ^ ^ 

for k = 0,2, . . . , D, and is defined in Eq. (j3.14p . For the off-diagonal component we have 

\ o/b-^ 2Dvr^/2+ia[a(a-r)/r/]^' ^' ^ 

The leading terms do not depend on the mass and have opposite signs for Dirichlet and non- 
Dirichlet boundary conditions. In particular, for a minimally coupled Dirichlet scalar field 
the energy density near the sphere is negative and the energy flux is directed away from the 
boundary. 

In Fig. [H for D = 3 dS space, we have plotted the boundary-induced part in the VEV 
of the energy density for a scalar field with the Dirichlet boundary condition as a function of 
the proper distance form the center of the sphere measured in units of the dS curvature scale 
a. (In figures below we plot the VEVs for both interior and exterior regions. The calculations 
for the exterior region are presented in Sect. El) For the sphere's proper radius (in units of a) 
we have taken the value a/r/ = 1. The left /right panel corresponds to comformally /minimally 
coupled scalar fields. The numbers near the curves are the values of the parameter ma. These 
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values are chosen in a way to have cases of both real and purely imaginary u. For a minimally 
coupled field the dependence of the VEVs on the mass is weak for points near the sphere. This 
follows from the asymptotic expansion (|4.2U|) . For a conformally coupled field the leading term 
vanishes and the subleading terms in the asymptotic expansions depend on the mass. Note 
that for conformally and minimally coupled fields the graphs have different scales. For non- 
Dirichlet boundary conditions, the VEV of the energy density near the boundary has opposite 
sign compared with the Dirichlet case. 




1 2 3 4 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

rjt] r/r] 



Figure 1: The boundary- induced part in the VEV of the energy density as a function of the 
proper distance from the center of the sphere for conformally (left plot) and minimally (right 
plot) coupled D = 3 scalar fields with the Dirichlet boundary condition. The numbers near the 
curves are the corresponding values of the parameter ma and for the sphere's proper radius we 
have taken a/rj = 1. 

Fig. [2] displays the dependence of the boundary-induced part in the energy density, for a 
fixed proper distance from the center of the sphere with the Dirichlet boundary condition, as 
a function of the mass measured in units of the dS energy scale. As in Fig. [H the left and 
right panels are for D = 3 conformally and minimally coupled fields respectively and the sphere 
radius is taken a/r] = 1. The numbers near the curves are the values of the ratio r/r/. For large 
values of ma the VEVs exhibit an oscillatory behavior. 

Fig. [3] shows the energy flux as a function of the proper distance from the center of the sphere 
with the radius a/rj = 1 for D = 3 conformally and minimally coupled fields with the Dirichlet 
boundary condition. Again, for a minimally coupled scalar field the dependence on the mass is 
weak, which directly follows from the asymptotic expression (j4.2ip . For both cases of minimally 
and conformally coupled fields the energy flows out of the sphere. For non-Dirichlet boundary 
conditions the energy fiux has the opposite sign. 

The dependence of the energy fiux on the mass of the field is depicted in Fig. [5] for a fixed 
value of the proper distance from the sphere's center. As before, the case of a D = 3 scalar field 
with the Dirichlet boundary condition is considered. For a conformally coupled massless field 
the energy flux vanishes 

In this section, we have investigated the VEV of the bulk energy- momentum tensor inside 
a spherical boundary. For a scalar flelds with Robin boundary conditions, it was shown in Ref. 
|49j that in the discussion of the relation between the mode sum energy, evaluated as the sum 
of the zero-point energies for each normal mode of frequency, and the volume integral of the 
renormalized energy density it is necessary to include in the energy a surface term concentrated 



14 





Figure 2: The boundary-induced part in the VEV of the energy density as a function of the field 

mass for a fixed proper distance from the center of the sphere in Z? = 3 dS spacctime. The left 
and right panels correspond to conformally and minimally coupled fields. The numbers near the 
curves are the values of the proper distance from the sphere's center in units of a. 
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Figure 3: The energy flux as a function of the proper distance from the center of the sphere for 
D = 3 conformally and minimally coupled scalar fields (left and right plots respectively) with 
the Dirichlet boundary condition. The numbers near the curves correspond to the values of the 
parameter ma. 
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Figure 4: The energy flux as a function of the mass for a fixed proper distance from the center of 
the sphere with the proper radius a/r] = 1. The left and right panels are for D = ?> conformally 
and minimally coupled fields with the Dirichlet boundary condition. The numbers near the 
curves are the values of the ratio r/r] (proper distance from the sphere's center in units of a). 

on the boundary (see also Ref. [50]). An expression for the surface energy-momentum tensor 
for a scalar field with a general curvature coupling parameter in the general case of bulk and 
boundary geometries is derived in Ref. [51j . By making use of this expression and the mode 
functions from Eq. (12. Qh . we can evaluate the mode sum for the surface energy- momentum tensor. 
The latter is located on the sphere and, in addition to the divergences of the boundary-free dS 
spacetime, contains surface divergences. The corresponding renormalization procedure is similar 
to that for the total Casimir energy and will be discussed elsewhere. 

5 Wightman function in the exterior region 

Now we turn to the evaluation of the VEVs in the region outside a spherical shell with radius a. 
The boundary condition for a scalar field is in the form (j2.4p . As we did for the interior region, 
first we consider the Wightman function. This function is evaluated by using the mode sum 
formula ()2.5p . In the exterior region, r > a, the radial parts of the mode functions are given by 
Eq. (|2.7p , where the ratio of the coefficients in the linear combination of the Bessel and Neumann 
functions is determined from the boundary condition ()2.4p . Note that the boundary condition 
can be written in the covariant form (1 + I3n^\/i)(p = 0, with being the outward normal (with 
respect to the region under consideration) to the boundary normalized as gun^n^ = — 1. For the 
interior and exterior regions we have n' = (r//a)(5[ and = — (ry/a)(5[, respectively. Comparing 
with the boundary condition in the form (12. 4p . we see that B/A = (r]/a)f3 for the interior region 
and B/A = —{r]/a)l3 for the exterior one. Hence, when we impose Robin boundary condition 
with the same coefficient /3 then the ratio B/A has opposite sign for the exterior and interior 
region. In general, (3 could be different for these regions. 

In the exterior region, for the mode functions realizing the Bunch-Davies vacuum state, one 
finds 

^„{x) = C,-^Lj^Hi'){Xrj)g^{Xa,Xr)Y{mk;'&,4>), (5.1) 

where we have defined 

g^iXa, Xr) = Y^{Xa)J^{Xr) - J^{Xa)Y^{Xr), (5.2) 
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with ^ A < cxD and with the barred notation (|2.2U|) . The normahzation coefficient is determined 
from the condition ()2.12p . where now the integration goes over the exterior region and in the 
right-hand side the delta symbol for A is understood as the Dirac delta function 6{X — A'). 
In evaluating the normalization integral over the radial coordinate we note that this integral 
is divergent for A = A' and, hence, the dominant contribution comes from large values of r. 
By making use of the asymptotic formulas for the Bessel and Neumann functions for large 
arguments, from Eq. (12.12P one gets 

where we have used Eq. (j2.14p . 

Substituting the functions (15. ip into the mode sum (I2.5p , for the Wightman function we find 
the expression: 

^ 1=0 

°°dAA «'-'^°'^;>f'^<^°-Y) gW(A,)gg'(A,-). (5.4) 

Similar to the case for the interior region, this function can be written in the decomposed form 
()2.22p . In order to extract the boundary- induced part explicitly, we subtract from Eq. (j5.4p the 
function W(is{x,x'), written in the form (|2.23p . For the further evaluation of the difference we 
use the relation 

g^(Aa,Ar)g^(Aa,A/) . - -i V Ai^n'-^) ( Xr)H^J^ ( Xr') (5 5) 

As a result, the boundary- induced part is presented in the form 



oo 



X A Jii^H<'>(A,-)Jft'>(Ar')JfJ"(A^)^?i?(A,'). (5.6) 

(Aa) 

Assuming that the function Hj^\z), {hI;^\z)) has no zeros for < argz ^ ■k/2 (—-k/2 ^ argz < 
0), we rotate the integration contour by the angle tt/2 for the term with j = 1 and by the angle 
— 7r/2 for j = 2. Introducing the modified Bessel functions, for the boundary-induced part of 
the Wightman function in the exterior region we get the expression 

(„„>-)D/2 ^ poo 

dx X 



W,{x,x') = , ^ Y.^2l + n)C-'\co.e) / 



#r4^M(^2;)K^(r'x) [I,{xr,')K,{x1^) + I^,{xr^)K,{xr,')\ . (5.7) 

Comparing with Eq. (j2.24p . we see that the expressions for the boundary induced parts in the 
exterior and interior regions are related by the interchange ^ I^. 
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6 Vacuum expectation values in the exterior region 



6.1 Field squared 

First we consider the VEV for the field squared. It is presented in the decomposed form ()3.ip . 
with the boundary-induced part given by 

= - of/V2 EA r dxx'-^MpIl!lKlirx/rj)F^ix), (6.1) 
TrSD{r/vr ^ Jo K^{ax/r]) 

where Fu{x) is defined by Eq. (14. Sp . The time dependence in this expression appears through 
the proper radius of the sphere and the proper distance from the center. For a conformahy 
coupled massless field one has the relation (13. 6p . where the expression for ((/?^)m is obtained 
from Eq. (j3.7p by the interchange ^ I^. For points near the sphere the leading term in 
the asymptotic expansion over the distance from the sphere is in the form (I3.13P with (a — r) 
replaced by (r — a). In particular, near the sphere the VEV has the same sign for the interior 
and exterior regions. 

Now we turn to the asymptotic behavior of the boundary-induced VEV at large distances 
from the sphere, assuming that r/rj ^ 1 for fixed a/rj. In this limit, the dominant contribution to 
the integral in Eq. ()6.ip comes from the region near the lower limit of the integration, x < r//r, 
and we can use the relation 

for z ^ 1 . By using also Eq. (13. 9p , the integral in Eq. (16. ip involving the square of the Macdonald 
function is evaluated with the help of a formula from Ref . [H] . The dominant contribution comes 
from the mode with / = 0. The behavior of the boundary-induced VEV at large distances is 
qualitatively different for positive and imaginary values of the parameter i^. In the first case, 
the leading term is given by the expression 



2a^-^ r(z^)r(n/4 + 1 - u) 
ttSd nr2(n/2)r(l - zv) 

A + nB/2 (a/7?)^-2 
A-nB/2 {r/r]f 



xF (3n/4 + 1 - .) .„\"(^_,,„„ (6.3) 



and the VEV is a monotonically decreasing function of the radial coordinate. 

At large distances from the sphere and for imaginary v, v = i\v\, the VEV behaves as 



Co{u)a^-^ A + nB/2 
T,D/2+iY^D/2- I) A-nB/2 
(a/n)^-"^ 

X lr/ii)2D-2 Hr/V) + M, (6.4) 

where the coefficient Co and the phase <j)o are defined by the relation 

Co{i^)e''^° = f^^r (3n/4 -iy + l)r (n/4 - u + I) . (6.5) 
In particular, for D = 3 one has 

,2 1 



13/16 -z^ + z.^ 

Co{u) = 7rV2 -n: . 

|z^|coshi/2(2^|z.|) 



3.6) 
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As is seen, for imaginary v the behavior of the boundary-induced VEV is damped oscillatory. 
For a fixed value of r, the VEV behaves as (v9^)b ~ e~^*/° cos(2|i/|t/a + (/)q), where (/)q = 

ln(r/a) + (/)o. The oscillation frequency is increasing with increasing mass of the field. 

It is of interest to compare the behavior of the VEV at large distances to the corresponding 
behavior for a spherical boundary in Minkowski spacetime. In the latter case, for a massless field 
the decay is of power-law: (v'^)m ~ a^~^/r^^~^. For a massive field, assuming that mr S> 1, 
the VEV behaves as 

(9^ )m - - , e ^n_i } Di-^ -, (6.7) 



4S/)r^~i ^/mr ^-^ K^{am) 

and it is exponentially suppressed. In the case of dS spacetime, for a fixed value of ma < 1, the 
decay at large distances is of power law (monotonic or oscillatory) for both massless and massive 
fields. So, we conclude that the curvature of the background spacetime essentially changes the 
behavior of the VEVs at distances larger than the curvature radius of the background spacetime. 

6.2 VEV of the energy- momentum tensor 

Now we turn to the VEV of the energy-momentum tensor outside the sphere. This VEV is 
presented in the decomposed form (|4.ip . By using the formula (j4.3p and Eq. (jG.ip for the VEV 
of the field squared, the diagonal components of the boundary-induced part are presented in the 
form (no summation over i) 



2tt bDa^+^ '7^ Jo K^{ax) 



1=0 

^2 



X {G'^K^[rx)]F,{xri) + 2Kl(rx)Fl{xr,)] , (6.8) 

where the functions F?(?/) and G\ [/(y)] are defined by Eqs. (|4.6p and (j4.7p . For the off-diagonal 
component we have the expression 



X [2rxKl{rx) - {D - 2)K^{rx)\ [{i - 1/4) r]d^ + i] F^{xr]). (6.9) 



These components are written in the coordinate system {T,r,^),(p). They are related to the 
components in the system (t, r, i?, (/>) by Eq. (j4.10p . 

The vacuum energy in the region r ^ rg > a is given by the expression (I4.14p with the 
integration over this region. Making use of the expression for the energy density from Eq. (j6.8p . 
one finds the formula 

^(b) _ rdxx'-^I^^^^^ 

D-2 



X 



{{l-4()F,ixri/ro)K^{x) 



xK'^(x)-^^K^{x) 



2 



-F^ixv/ro) [x^K';^{x) - {x^ + fi^)Kf^{x)] }. (6.10) 

This energy is related to the energy flux and to the vacuum stresses in the exterior region by 
Eq. (j4.15p . Similar to the interior region, the integrated energy diverges in the limit rg — )• a. 

Let us consider the asymptotics for the VEV of the energy-momentum tensor. For points 
near the sphere, in a way similar to that for the interior region we find for r — )• a (no summation 
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over k) 



DTiiD + l)/2)iC-CD)KB 
2D Tr(D+i) /2 [a{r -a)/ 7]]^+^' 



for k = 0,2, . . . , D, and 



iri\ ~ {D-i)ri{D + i)/2){C-CD)KB 

^ 1)^ - 2DTr{D+i)/2a(r-a)D{a/r])D+^ ' ^ ' 

/Ti\ ^, Dr{D/2){C-CD)KB 
^^o^b-^ 2^7r^/2+i«[a(r - a)/r?]^ • ^ ' 

As we see, for a scalar field with non-conformal coupling the components {T^)h, k = 0,2, . . . , D, 
have the same sign near the sphere for exterior and interior regions, whereas the radial stress 
and the off-diagonal component change signs. 

Now we consider the limit of large distances from the sphere, assuming that r/rj ^ 1 for 
fixed a/r]. Similar to the case of the field squared, discussed in the previous subsection, for 
positive values of v to the leading order we get: 

^ 27r^/2+irp/2 - l)a^+i A + {D - 2) B/a ^ ' 



In this expression, 

and (no summation over A;) 



DV{v) 



,^D-2v 
{D + l-2v)i 



r(i- 

xr(3n/4 + 1 - i/)r(n/4 + 1 - z^) , (6.14) 



Bl{y) = ^B^oi'^), k = l,2,...,D, 

Bl{v) = -2 ^~^~^ i?oV)- (6.15) 

In this case the boundary-induced VEV decays monotonically with increasing proper distance 
from the sphere. At large distances the vacuum stresses are isotropic. If we denote the corre- 
sponding effective pressure by P\^, P\, = —{T^)]^, then the equation of state is of the barotropic 
type: = -{2v / D){T^)^. 

At large distances from the sphere and for imaginary v, the leading term in the asymptotic 
expansion is in the form 

, ^ (^/^)^+i-?? {a/r)D-^C^{u) A cos[2\y\ ln(r/r?) + 0f ] 
^Wb^ vr«/2+ir(Z)/2 - l)a^+i A + {D-2)B/a ' ^ ' 

where Cf{v) and the phases (j)^ are defined by the relation 

BKu) = Cf{u)e''^'. (6.17) 

As in the previous case, to the leading order the vacuum stresses are isotropic. The oscillations 
in the energy density and the stresses are shifted by the phase 7r/2. For a fixed comoving radial 
distance r, the VEVs oscillate like (7;^)b ~ e"(^+^"'^»*)*/" cos(2|zy|t/a + </>f ). For the Neumann 
boundary condition, the leading terms, given by Eqs. (j6.13p and (j6.16p . vanish and it is necessary 
to take the next terms in the asymptotic expansions. As a result the asymptotic expressions 
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for the Neumann boundary condition contain an additional factor (a/r)"^ compared with non- 
Neumann boundary conditions. In Figs. [IHll the boundary-induced part of the energy density 
and the energy flux in the exterior region are depicted as functions of the proper distance from 
the sphere center and of the mass, for D = 3 minimally and conformally coupled fields with 
the Dirichlet boundary condition and for the sphere's proper radius corresponding to a/r] = 1. 
Both cases of real and imaginary v are considered. The oscillatory behavior at large distances in 
the latter case is plotted separately in the inset. The first zero decreases with increasing mass, 
whereas the oscillation frequency increases. 

7 Further generalizations 

The approach described above may be used for the evaluation of the Casimir densities for more 
general background geometries, in particular, in models where the spherical boundary separates 
spacetime regions with different geometries. Here, as an example, we consider a model in which 
the spacetime is described by two distinct metric tensors in the regions outside and inside the 
spherical boundary with radius a. The line element in the exterior region will be taken in the 
form (|2.2p . We will assume that inside the sphere the spacetime geometry is regular and is 
described by the general spherically symmetric line element 

ds^ = aV-2[e2"('-)dt2 _ e2-Mdr2 - e^^-^df]!,], (7.1) 

where the functions u{r), v{r), w{r) are continuous at the core boundary: 

u{a) = v{a) = 0, w{a) = ln(a). (7.2) 

Here we assume that there is no surface energy-momentum tensor located at r = a and, hence, 
the derivatives of these functions are continuous as well. Note that by introducing the new radial 
coordinate r = e'^^^^ with the sphere's center at f = 0, the angular part of the line element ()7.ip 
is written in the standard Minkowskian form. With this coordinate, in general, we will obtain a 
non-standard angular part in the exterior line element ()2.2p . From the regularity of the interior 
geometry at the sphere center one has the conditions u{r),v{r) — t- 0, and w{r) ~ Inf for f — t- 0. 
The curvature scalar for the metric corresponding to (17. 1[) is presented in the form 

R = {a/7]y^ R + D{D + l)a-2e-2"«, (7.3) 

where 

R = _2e-2^ [-u" + n'2-nV + (L>-2)(L'-l)u;'V2 

+{D - 1) {w" + w'^ + w'u' - w'v')] +{D- 2){D - 1)6-^"" (7.4) 

is the curvature scalar for the static metric corresponding to the line element inside the square 
brackets in (j7.ip . Explicitly writing the field equation (j2.3p for the metric (j7.ip . it can be seen 
that the time and radial variables are separated in two special cases: for u(r) = or for a 
massless field. We will consider the first case. The discussion for the second case is similar. 

Taking u{r) = and assuming the Bunch-Davies vacuum state, the mode functions in the 
region inside the sphere, r < a, are written in the form 

V^ix) = ??^/2^«(A??)/Kr)y(mp;T9,</)), (7.5) 

where the radial function is a solution of the equation 

fl'ir) + [{D - l)w' - v']f[{r) + e^^ [A^ -^R-l(l + D- 2)e-^^] fi{r) = 0. (7.6) 



21 



and 

R = -2{D - l)e-2^ {w" + Dw'^/2 - w'v') + {D - 1){D - 2)e-^'^ . (7.7) 

The solution of the radial equation (j7.6p regular at the origin we will denote by i?/(r, A). 
Near the center this solution behaves like f'. Note that the parameter A enters in the radial 
equation in the form A'^. As a result the regular solution can be chosen in such a way that 
i?/(r, — A) = const • Ri{r, A). Now for the radial part of the eigenfunctions one has 

f ( \ = l for r < a 

\ [AiJ^{\r) + BiY^{\r)] for r > o " ^ ' 

The coefficients Ai and Bi are determined by the conditions of continuity of the radial function 
and its derivative at r = a: 

Ai = |a^/2-ii?,(a,A)i;(Aa), 

Bi = -^a''/^-'Ri{a,X)Jf,{Xa), (7.9) 



with the notation, for a cylinder function F{z), 

F(z) = zF'(z) 



D R',(a,z/a) 



2 Ri{a,z/a) 



F{z), (7.10) 



where i?;(a,A) = drRi{r^X)\r=a- Due to our choice of the function Ri{r,X), the logarithmic 
derivative in formula ()7.10p is an even function of z. Hence, in the region r > a the radial part 
of the mode functions becomes 

Mr) = |(a/r)^/2-ii?,(a,A)V(Aa,Ar), (7.11) 

with the notation 

h^iXa, Xr) = Yf,{Xa)Jf,{Xr) - Jf,{Xa)Y^{Xr). (7.12) 
The normalization condition (I2.12P is written in terms of the radial eigenfunctions as 

*/M/,(.-.A)/,(..A') = j|i^e*-)-/=. (7.13) 

where ro is the value of the radial coordinate r corresponding to the sphere center and gr is the 
radial part of the determinant g. As the integral on the left is divergent for A' = A, the main 
contribution in the coincidence limit comes from large values r. By using the expression ()7.1ip 
for the radial part in the region r > a and replacing the Bessel and Neumann functions by the 
leading terms of their asymptotic expansions for large values of the argument, it can be seen 
that from (I7.13P the following result is obtained: 

o2. ^^ a2-^ai-^Ae^('^--*W2 ^^^^^ 
Rf(nX) = ^ T- (7.14) 

By taking into account relation (I7.14p , the mode functions in the exterior region are presented 
in the form ^ 

'^-(^) = c.^;J7^^(^)(Ar?)V(Aa,Ar)y(mfc;i9,</>), (7.15) 

where 

12 



c: 



4a^-iiV(mfc) 



[jl{Xa) + Y^{Xa)]-\ (7.16) 
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Now comparing with (jS.ip . we see that the mode functions (|7.15|) are obtained from the corre- 
sponding functions for the sphere with Robin boundary condition by the replacement 

A D R',(a,X) 

As a result, the Wightman function in the exterior region is given by expression (j5.4p with this 
replacement. 

Further evaluation of the Wightman function and the Casimir densities is the same as in the 
case of the Robin sphere, described in Sections [5] and [H The Wightman function in the region 
r > a is presented in the decomposed form (j2.22p . where now the part induced by the geometry 
()7.ip with u[r) = in the region r < a is given by (j5.7p with 

The factor e'^*/^ in the argument of the function Ri is related to the complex rotation we have 
used after formula ()5.6p . Similarly, the expressions for corresponding parts in the VEVs of the 
field squared and the energy-momentum tensor are presented as (j6.ip and (j6.8p with the ratio 
of the coefficients from (jT.lSp . 



8 Conclusion 

In this paper, we have investigated the Casimir densities induced by a spherical boundary in dS 
spacetime for a massive scalar field with general curvature coupling parameter. On the sphere 
the field obeys Robin boundary condition with coefficients, in general, different for the interior 
and exterior regions. We have assumed that the field is prepared in the Bunch-Davies vacuum 
state which is dS-invariant in the boundary-free spacetime. In free field theories all properties of 
a quantum field are encoded in two-point functions. We have computed the Wightman function 
and the Green's function both in the region interior to the sphere and in its exterior. The 
VEVs for the field squared and the energy-momentum tensor are obtained from these two-point 
functions in the coincidence limit. In addition, the Wightman function determines the response 
of particle detectors of the Unruh-DeWitt type. In a similar way other two-point functions can 
be investigated. 

For the evaluation of the Wightman function we have employed the mode summation method. 
In the region inside the spherical boundary, the mode functions for a scalar field, realizing the 
Bunch-Davies vacuum state, are given by Eq. ()2.9p . The eigenvalues for A are quantized by the 
boundary condition on the sphere and they are the solutions of Eq. ()2.10p . The corresponding 
mode sum for the Wightman function is given by Eq. ()2.17p and contains the summation over 
these eigenvalues. The latter are given implicitly and Eq. (j2.17p is not convenient for the further 
evaluation of the VEVs in the coincidence limit. In order to obtain a more workable form, we 
have used the generalized Abel-Plana formula for the summation over the eigenvalues of A. This 
allowed us (i) to extract explicitly the boundary-free Wightman function and (ii) to present 
the sphere-induced part in terms of an integral rapidly convergent in the coincidence limit (for 
points away from the boundary). The local geometry away from the sphere is the same as in 
the boundary-free dS spacetime and the renormalization for the VEVs of the field squared and 
the energy-momentum tensor is reduced to that for the boundary-free dS spacetime. The latter 
is well investigated in the literature. In the same way we can evaluate the Wightman function 
in a general vacuum state corresponding to the mode functions (|2.7p with a linear combination 
of Hankel functions. The corresponding expression for the boundary-induced part has the form 
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similar to Eq. (|2.24|) with a more general bilinear combination of the modified Bessel functions 
of order v. 

By using the representation of the Wightman function, the VEVs of the field squared and 
the energy-momentum tensor are decomposed into the boundary-free and boundary-induced 
parts. For the region inside the sphere the latter are given by Eqs. (|3.3p . (j4.4p . and (j4.9p . for the 
field squared and the energy-momentum tensor, respectively. An important feature is that the 
vacuum energy-momentum tensor has an off-diagonal component (14. 9p which describes energy 
flux along the radial direction. With dependence on the boundary condition and the mass of 
the field, this flux can be either positive or negative. The boundary induced VEVs for both 
field squared and the energy-momentum tensor depend on time through the proper radius of 
the sphere and the proper distance from the sphere's center. This property is a consequence of 
the maximal symmetry of the Bunch-Davies vacuum state. 

We have explicitly checked that the boundary-induced part in the VEV of the energy- 
momentum tensor obeys the trace relation (I4.11|) and the covariant conservation equation. In 
particular, the energy-momentum tensor is traceless for a conformally coupled massless field. 
The trace anomaly is present in the boundary-free part only. For a conformally coupled mass- 
less scalar field the flux vanishes. In this case we have simple relations, Eqs. (j3.6p and (j4.13p . 
between the boundary-induced VEVs for spherical boundaries in dS and Minkowski spacetimes. 
The latter are consequences of the conformal relation between the problems in Minkowski and 
dS spacetimes. Divergences are found in the VEVs as the surface of the sphere is approached. 
The leading terms in the corresponding asymptotic expansions in terms of the distance from the 
boundary are given by Eqs. (I3.13p . ()4.20p and ()4.2ip . Written in terms of the proper distance 
from the boundary, for a non-conformally coupled field these leading terms coincide with the 
corresponding terms for a spherical boundary in Minkowski spacetime. This is a consequence of 
the fact that for points near the boundary the dominant contribution to the VEVs comes from 
the modes with small wavelengths which are not influenced by the gravitational field. 

In the region outside a spherical shell, the eigenvalues for the quantum number A are contin- 
uous and the modes of the field realizing the Bunch-Davies vacuum state are given by Eq. (|5.ip . 
We have explicitly extracted from the Wightman function the part corresponding to dS space- 
time without boundaries. The boundary-induced part is given by Eq. (j5.7p for the Wightman 
function and by Eqs. (16. ip and (16. Sp . ()6.9p for the field squared and the energy-momentum ten- 
sor, respectively. General formulas are simplified in the asymptotic regions. For points near the 
boundary the leading terms in the expansions over the distance from the sphere are given by 
expressions ()6.1ip and (I6.12p . In this region the energy density and the azimuthal stresses have 
the same sign for the exterior and interior regions, whereas the radial stress and the energy flux 
have opposite signs. In particular, for a minimally coupled scalar field the energy flows away 
from the boundary for Dirichlet boundary condition and toward the boundary for non-Dirichlet 
boundary conditions. 

Most interesting is the behavior far from the sphere; qualitatively different behavior occurs 
depending on the sign of = /2 — £,D{D -|- 1) — rr?G? ^ where D is the number of spatial 
dimensions, ^ is the conformal parameter, m is the mass of the scalar field, and a is the curvature 
scale. When v is positive, the mean field squared and the energy-momentum tensor fall off as 
a power, while when v is imaginary, the large distance behavior is damped oscillatory with 
the amplitude decaying as (rj/r)"^^^'^^ for the field squared and the diagonal components of 
the energy-momentum tensor. For a scalar field with Neumann boundary condition the VEVs 
at large distances are suppressed by an additional factor (rj/r)'^ compared with the case of 
non-Neumann boundary conditions. Note that the behavior of the VEVs at distances larger 
than the curvatures scale of the background spacetime is completely different from the case of 
a spherical boundary in Minkowski spacetime. In the latter case and for a massive field the 
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boundary-induced VEVs at distances larger than the Compton wavelength of the scalar particle 
are exponentially suppressed by the factor e~'"(''~"). For the problem in dS spacetime, under 
the condition m < a~^, the decay of the VEVs is power-law. Exponential damping can occur 
only in an intermediate region r ^ 1/m. 

In Section[7]we have generalized the expressions of the Casimir densities in the exterior region 
for a class of spherically-symmetric metrics in the region r < a described by the line element 
()7.ip . The geometry in the exterior region is given by the dS line element (12. 2p . A special 
case with u{r) = is considered when the time and radial variables in the field equation are 
separated. The VEVs in the exterior region are decomposed as the sum of pure dS part and the 
part induced by the geometry in the interior region. We have shown that the expressions for the 
latter are obtained from the corresponding expressions outside the Robin sphere, investigated 
in Sections [5] and m taking the ratio of the Robin coefficients in the form ()7.18p . where Ri{r,X) 
is the solution of the interior radial equation regular at the origin. 

At the end, we would like to emphasize that the main subject of the present paper is the 
investigation of the local characteristics of the vacuum, the VEVs for the field squared and the 
energy-momentum tensor, at the points away from the boundaries. They do not contain surface 
divergences and are completely determined within the framework of standard renormalization 
procedure in quantum field theory without boundaries. We expect that similar results would be 
obtained in the model in which instead of externally imposed boundary condition the fluctuating 
field is coupled to a smooth background potential that implements the boundary condition in a 
certain limit [52]. 
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A Boundary-free part of the Wightman function 

For the further evaluation of the boundary-free part, given by Eq. ()2.23p . firstly we apply 
Gegenbauer's addition theorem for the Bessel functions (see, for instance, Ref. ^38j) to the series 
over I. This gives the following expression for the Wightman function: 

W (x x') - Tre^^'--'-*)"/^ 2-"/2-iai-^ (r?V)^/' 
~ r{n/2)nSD (r2 + r'2 - 2rr' cos Of^^ 

POO 

X / dzz''/'^+'^Hl^\zT])HlV{zr]')J,,/2izy^r^ +r'^ -2rr' COS 6). (A.l) 

As the next step, we write the product of the Hankel functions in terms of the Macdonald 
function: 

gi(.-.*W2^(i)(^^)^{2)(^^/) = -^K,{-izr])K,{izr]'), (A.2) 
and use the integral representation [H] 

K.iZ)K,A.) ^ i Z;;.,^"^.-".— — .e.p (_| _ ^) (A,3) 
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for the product of the Macdonald functions. Substituting Eqs. ()A.2p and ()A.3p into Eq. (jA.ip . 
we first integrate over x and then with respect to w, with the result 



where 

n(x, x') = 1 + ^^^^ —, . (A.5) 

In deriving Eq. ()A.4p we have assumed that |m| < 1. The integral in Eq. ()A.4p is expressed 
in terms of the associated Legendre function x')) (see Ref. [53]). Expressing this 

function through the hypergeometric function, after some transformations, we get the final 
expression for the Wightman function in dS spacetime (for two-point functions in boundary-free 
dS spacetime see Ref. 



a'-^ TiD/2 + ,^)TiD/2 

Wds{x,x) - 



(47r)(^+i)/2 r((i? + l)/2) 

^ /D D D + l l + u(x,x')\ , 
x.F,^- + .,--.;^; i^j. (A.6) 

Note that, if we denote by X{x) the coordinates in the higher-dimensional embedding space 
for dS spacetime, then one can write u{x,x') = 1 + [-'^(x) — X{x')]'^ /{2a'^). The property that 
the Wightman function depends on spacetime points through [X(x) — X{x')]'^ is related to the 
maximal symmetry of the Bunch-Davies vacuum state. 



B Minkowski spacetime limit 

In this appendix, for the boundary-induced part in the Wightman function, we explicitly demon- 
strate the limiting transition to the geometry of a spherical boundary in Minkowski spacetime. 
In the Minkowski spacetime limit one has a — )• 00 and the modulus of the order of the modified 
Bessel functions in Eq. ()2.22p is large, u ~ ia, a = ma S> 1. In addition, we have rj ^ a — t. We 
make use of the uniform asymptotic expansions for the modified Bessel functions for imaginary 
values of the order with large modulus. For z < 1, the leading terms in these expansions have 
the form (see, for example, Ref. |54] ) 



^e-''-/\os[af{z)-7T/A], 



2g<T7r/2 

h^iaz) + I-^Aaz) ~ --^sm[af{z) - 7r/4], (B.l) 



where 



m = In ^^^^-^ - yr^. (b.2) 



z 

In the case z > 1 one has the asymptotics 

— g-o-7r/2-(7g(z) 



Kia{(Tz) 



2a (z2_ 1)1/4 



2 ^an/2+ag{z) 

liaiaz) + I-^a{cTz) 7^7^ 77777, (B.3) 
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with ^ 

g{z) = — arcsecz + \/z^—l, g'{z) = —\l — 1. (B.4) 

From Eqs. (|B.ip and ()B.3P it follows that the dominant contribution to the boundary-induced 
part of the Wightman function in Eq. (j2.24p comes from the integration range x > m. In this 
range we have 

, X T / X / /X coshima\q(z) — q(z')]} , 
h{xr^')K,{x^) + I.,{x7^)K,{xv') « , I AimV 2 iU/4 ' (I^-^) 

where z = xrj/ma and z' = xr]' /ma. By taking into account that r//a ~ 1 — t/a and r/'/o ~ 
1 — f'/a, it can be seen that 



g{z)-g{z')^^{x/mf -l{t' -t) la. (B.6) 
Now, combining this with Eq. (jB.Sp . we get 



T , i\r^ , \ T / , i\ cosh(At\/x2 - m^) 

ay x^ — 

with = t' — t. Substituting this into the expression for the boundary-induced part of the 
Wightman function, Eq. (I2.24p . to the leading order one finds 

Wb(x,x') « -^^^^ — y(2^ + n)C7"/2(cos0) / dxx 

imbD j-i Jm 



1=0 



Knixa) coshfAt-v/x^ — m?) ^ , , -r , , 
X / \ '-lJxr)lJxr'). (B.^ 



The expression in the right-hand side coincides with the Wightman function inside a spherical 
boundary in the Minkowski bulk [26j. The Minkowski space limit for the Wightman function in 
the exterior region is considered in a similar way. 



C Green's function 



In the text, we computed the Wightman function in order to calculate the VEVs of the field- 
squared and the energy-momentum tensor. Of course, these quantities can equally well be 
computed from the Green's function for this problem. The calculation of the latter is similar to 
that given for the Wightman function, but there are some points of interest, so we sketch the 
derivation here. 

Because the nontrivial structure in the dS geometry lies in the time dependence, it is natural 
to isolate that by constructing a reduced Green's function in the variable ?y, ??' {n = D — 2): 

= ,, ^(2/ + n)C;/^(cosg)^-r^(Aa)J^(Ar)J^(ArO/fa7?0. (C.l) 

I — U A 

Here the notation is as in Sec. [21 and in particular, the eigenvalues A are the roots of Eq. (j2.10|) . 
The reduced Green's functions satisfies 

f{v,v') = -{lY''s{v-v'). (C.2) 



1-D d 
7] drj 



+ 



am 



+ CD{D-1) 
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The solution of this equation corresponding to the Bunch-Davies vacuum (the analog of outgoing 
wave solutions) is 



(C.3) 



which uses the Wronskian (j2.14p and ?]> (r/<) is the greater (lesser) of rj, r]'. 

To resolve the difficulties with oscillatory integrals and an implicit equation for the eigen- 
values, we can again use the generalized Abel-Plana formula (j2.19p . which leads to the breakup 
of the Green's function into a free part referring only to the dS background, and a term which 
exhibits the effect of the sphere: 



where 



Gds{x,x') 



G{x, x') = Gds{x, x') + Gb(x, x' 



(C.4) 



vr e 



4i q^-^hSd {rr'Y/"^ 



^(2Z + n)C"/^ (cos 



X / (iAAJ^(Ar)J^(Ar')F«(A^r>)F(?)(Ar<), 



(C.5) 



and 



Gb(x,x') 



1 a 



l-D 



(w 



^(2/ + n)Cf/^ (cos 



■ni nSu (rr')"/^ 



1=0 



I^{Xr)I^{Xr') [i^,(Ar?>)/_,(A?7<) + K,(Ar?<)4(Ar/>)](C.6) 



There is, in fact, no discontinuity associated with the last factor, since the quantity in square 
brackets here is 

[/_,(Ar/)/_,(Ar/') - /,(Ar?)/,(A??')] ■ (C.7) 

smvrz^ 

Of course, Gb(x,x') is simply iWb(a;,x') given in Eq. ()2.24p . 

The calculation of the Green's function for the exterior region proceeds similarly. Again, if 
we impose the boundary condition at the smallest value of r in the region, here at r = a, we 
have the radial function given in Eq. (j5.2p . which have continuum normalization 



1 

drrg^iXr)g^{Xr') = -<^(A - A')[J^(Aa) + Y^iXa)]. 



(C.8) 



Thus the exterior Green's function is 

in e* 



G{x,x') 



4 a^-^nSo 



^(2/ + L>-2)C"/^(cos0) 

F(i)(Ar?>)/72)(A7?< 



(C.9) 



/o J^^iXa) + Y^iXay 

This evidently leads to the boundary-dependent part given by Eq. (j5.7p multiplied by i. 

Incidentally, we might record here the imaginary rotations of the Hankel functions, since 
these are not given in the standard tables: 



Hi^\e-^''^x) 



e'^^^/^I^^ix) + -e-'-^'^/^K^ix 
vr 



(C.IO) 



for real or imaginary v, which terms are related by evident complex conjugation. 
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